Superfluidity of Bosons in Kagome Lattices with Frustration 
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In this letter we consider spinless bosons in Kagome lattice with nearest-neighbor hopping and 
on-site interaction, and the sign of hopping is inverted by insetting a n flux in each triangle of 
Kagome lattice so that the lowest single particle band is perfectly flat. We show that in the high 
density limit, despite of the infinite degenerate of the single particle ground state, interaction will 
select out the Bloch state at the K point of Brillouin zone for boson condensation at the lowest 
temperature. As temperature increases, we show that the single boson superfluid order can be easily 
destroyed, while an exotic triple-boson paired superfluid order will remain. This trion superfluid 
exist in a broad temperature window until the temperature is increased to the same order of hopping 
and then the system turns into normal phases. Finally we show that time-of-flight measurement of 
momentum distribution can be used to distinguish these three phases. 



Flat band models have attracted considerable theoret- 
ical interests recently [THIS], because the single particle 
energies are degenerate inside the flat band, therefore 
interactions play a dominative role in the many-body 
system and lead to many interesting quantum phases [TJ- 
13 . Among many different physical realizations of lattice 
models with flat band, Kagome lattices with only near- 
est neighbor hopping is perhaps one of the simplest. Re- 
cently, such a model has been realized experimentally us- 
ing optical lattices by the Berkeley group [T3]. However, 
the flat band in a normal Kagome lattice is the highest 
band. By fast shaking the optical lattices, one can invert 
the sign of hopping, which has also been demonstrated 
experimentally for triangular optical lattices [15 . This 
technique can be applied straightforwardly to Kagome 
lattices. When the sign of hopping is inverted, it is equiv- 
alent to insetting a it flux in each triangle of the Kagome 
lattice, and the flat band becomes the lowest band, as 
shown in Fig.[T} Thus, the kinetic energy is frustrated on 
such a Kagome lattice. 

In this letter we consider spinless bosons with on- 
site interaction on the Kagome lattice with frustra- 
tion, and we focus on the high density limit, where 
the system is away from strongly correlated phases like 
Mott insulator [3 15], Wigner crystal [5] and quantum Hall 
state[ini [TTJ. This high density limit corresponds to a 
real system with a Kagome optical lattice in the xy plane 
and weak confinement potential along the z direction. 
Therefore, each site is in fact a tube inside which bosons 
form a quasi-condensate, and tunneling couples different 
tubes that leads to a two-dimensional Josephson array 
described by a two-dimensional boson Hubbard model. 

Here we shall discuss whether and how bosons can con- 
dense in the hatband, and what exotic type of superfluid 
occurs at finite temperature. The main findings are 

(1) Though the lowest band of the single particle 
state is infinitely degenerate, interaction will select out 
a unique single particle state for bosons to condense at 
the lowest temperature. We show that this unique single 
particle state is the Bloch state at the K point. 



(2) As temperature increases, the single boson super- 
fluid order can be easily destroyed. However, a triple- 
boson paired superfluid order will remain until the tem- 
perature is of the order of hopping energy and the system 
turns into normal. Thus we predict a large temperature 
window in which the system has an exotic "trion super- 
fluid" phase. 

(3) We show that a simple time-of-fight (TOF) detec- 
tion of momentum distribution can distinguish K point 
conventional condensate, trion condensate and normal 
phase. 

It has been a long term efforts to search exotic bosonic 
state that has no single-boson superfluid order. In pre- 
vious studies, two-boson paired superfluids (or charge-4e 
superconductor as paired Cooper pairs) have been pro- 
posed in various systems, such as stripe superconductor 
[32], FFLO superconductor [TgH20]. spin-1 bosons [UJH] 
and bosons with Rashba spin-orbit coupling [25]. These 
boson paired superfluid phase emerges from melting of 
stripe order or internal spin order. But none of them have 
been observed yet. As far as we know, it is the first time 
that a triple-boson paired superfluid is proposed, and the 
underlying physical mechanics, that is, the kinetic en- 
ergy frustration due to the particular Kagome geometry, 
is also different from that of previous examples. These 
results may also be generalized to interacting bosons in 
other flat band models with geometric frustration. 

Band Structure and Mean-Field. The boson Hubbard 
model on the Kagome lattice with tt flux in each trian- 
gle is given by H = H t + Hjj, with the hopping term 
H t = tY2(ij\ tyPi + h.c., and the interaction term Hjj — 

U Yji ni{hi -l)-fiJ2i ™i> w here n t = b\bi. Here t is pos- 
itive [H] . In the momentum space, the hopping Hamilto- 
nian reads H t = J2k b\h{k)b k with b k = (b k A, bkB, bkc) 1 
and 

(0 cosfc3 cosk 2 \ 
cosfc 3 cosfci , (1) 
COS/C2 cosfci J 

where fc$ = k ■ a i: and a x = (1,0), a 2 = (— §> 2 )> 
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FIG. 1: (a) Kagome lattice, partitioned into A, B, C sublat- 
tices. (b) Band structure with inverted hopping sign, where 
r = (0,0), K = (27r/3,0) and M = (tt/2, 7r/2v^3). (c) Mean 
field energy landscape for different Up/t. Energy shifted by 
Eq = —2t + Up/3 (p = the boson density per unit cell). 




FIG. 2: (Color online.) (a-b) Phase configuration of the con- 
densate at r point (a) and K point (b), (c-d) random 3-color 
arrangement. Phases are denoted by colors: red = 1, green 
_ gZWs^ bme _ e 4W3 mar k out the vor ticity around 

each triangle. 



a 3 = ( — |) — 2 ) as shown m Fig-|l[a). The band struc- 
ture contains three bands: two are dispersive, and one is 
flat. The flat band is the lowest band, with a quadratic 
touching point between the lowest and the first excited 
band at T point, as shown in the Fig.jljb). 

Because of the infinite degeneracy of the single particle 
ground state, bosons can not condense if there is no inter- 
action. However, once the interaction energy is included, 
the single particle degeneracy will be lifted. Here we shall 
first discuss how the interaction select out a unique single 
particle state for boson condensation. At the mean-field 
level, we first consider the single particle state with trans- 
lational symmetry labelled by momentum k. Taking the 
mean-field ansatz — z = (za, z-q, zq) 1 with the nor- 
malization condition z^ z = p, where p is the boson filling 
per unit cell, we can minimize the whole mean-field en- 
ergy function E k [z] = z^h(k)z + £/(|z A | 4 + \z B \ 4 + |z c | 4 ) 
with respect to z, and denote the optimal mean field 
energy as E^ = min 2 E^ [z] . This mean field energy land- 
scape is shown in Fig.[ljc). 

Fig.[ljc) shows that the single particle degeneracy is 
lifted by the Hartree energy. This is because inside the 
flat band, only the single particle wave function of the 
K point has equal amplitude in three sublattices, and 
therefore only at K point the density is uniform and the 
repulsive interaction energy is the smallest. While nearby 
r point, a uniform density state can be achieved by form- 
ing a superposition between the lowest and the first ex- 
cited band, and such a superposition costs no energy at 
r point because of the quadratic band touching. Thus, 
the mean-field energy select out the T and K points. 

The Bloch wave function for T and K points deter- 
mined by mean-field energy read 

Mri) = ^(l,e ±2 "/ 3 ,e T2 "/ 3 ) T , 

, (2) 

?Mn) = 4,e ifc ^(i,-i,-i)T, 



with kx = (27r/3,0). In the real space, both two wave 
functions satisfy two conditions: (i) their densities are 
uniform, which minimize the mean-field interaction en- 
ergy and (ii) their phases follow "3-color arrangement" . 
Here "3-color arrangement" means that each pair of two 
neighboring sites take different phases out of 1, e 27 ™/ 3 and 
e -27™/3^ as depicted in Fig. [2] It is easy to show that once 
the "3-color arrangement" is satisfied, the kinetic energy 
is also minimized. "3-color arrangement" requires that 
three sublattices in each triangle must all take different 
phases, and one can define "vorticity" for each triangle. 
Vorticity +(— ) means that the phase increases clockwise 
(anti-clockwise). As shown in Fig.[2]ja) and (b), ipr is 
a "vorticity ferromagnetic" state and ipK is a "vorticity 
antiferromagnetic" state. 

In fact, once above condition (i) and (ii) are satisfied, 
both kinetic and interaction energy are simultaneously 
minimized at the mean-field level. The translation in- 
variance is not a necessary condition and can be released. 
In fact, one can find extensive number of states without 
translation invariance that satisfies both (i) and (ii), and 
two of such examples are depicted in Fig. ^ c-d). These 
states are locally identical to either K state or T state, 
and the domain wall costs no energy at mean-field level. 

Quantum Fluctuation and Order from Disorder. These 
extensive number of degenerate mean-field states can be 
further lifted by quantum fluctuations, because the zero- 
point energy of Bogoliubov phonons is different for differ- 
ent mean-field state. For a given mean-field configuration 
{hi) = Zi, its Bogoliubov Hamiltonian is 

H[zi] ^^(blk+h.c)-^^ 

1 (3) 
+ uY / (^ l \ 2 Hk + ^brbl + h.c.), 
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FIG. 3: (a) Bogoliubov spectra for T point (left) and K point 
(right) condensates at Up/t = 0.3. Momentum points are 
defined as A = (vr/3,0) and B = (tt/3, tv/S^S). (b) Distri- 
bution of the zero point energies of all degenerate mean-field 
configurations. The horizontal axes are zero-point energy A of 
degenerate mean-field state, and two dashed lines denote the 
zero-point energy for K and V point condensate. The vertical 
axes is the probability for having this zero-point energy, (c) 
J/t v.s. Up/t plot. 



where p — —2t + 2Up/3. Diagonalization of H[zi] leads 
to H[zi\ = I] m w m [zi] (jjnjm + 1/2), with 7 ro being Bo- 
goliubov boson operator. Thus the vacuum energy asso- 
ciated with this given condensate configuration equals to 
A N = | Em u m [zi] by setting ^ m % n = 0. 

For condensates at T and K points, the Bogoliubov 
excitations have quantum number k and they have well 
defined dispersion as shown in Fig.[3^a). One finds that 
the sound velocity c of K point condensate is smaller that 
that of r point condensate. This is also evidenced from 
the mean field energy landscape as shown in Fig.[ljc), 
where energy landscape changes less rapidly nearby K 
point than that nearby T point, indicating softer Gold- 
stone mode and more quantum fluctuations. Roughly 
speaking, the zero point energy is proportional to the 
sound velocity squared, i.e. A ~ c 2 , the degeneracy be- 
tween r point and K point configurations will thus be 
lifted by their different zero point energies. 

For a generic mean-field state in the degenerate man- 
ifold, the Bogoliubov spectrum has no well defined dis- 
persion because of the absence of translation symmetry. 
Here we randomly sample over 4000 configurations on a 
240-site Kagome lattice which have uniform density and 
satisfy "3-color arrangement" , and then calculate their 
zero-point energy numerically. We find that their zero- 
point energies all rest between T point condensate and 
K point condensate, i.e. V{zJ : h[K] < A[z 4 ] < A[r], 
as seen from Fig.[3]jb). Now the degeneracy has been 
completely removed through the order-by-disorder mech- 
anism. At sufficient low temperature, bosons will con- 
dense to K point (or its symmetry related points). 

Thermal Fluctuation and Phase Diagram. We have 
mentioned that all mean-field degenerate states can be 



locally viewed as either K-po'mt or T point, it is natu- 
ral to introduce an energy scale J = A[r] — A[K) as the 
typical energy cost for creating domain walls, which en- 
tirely comes from the fluctuation energy. As shown in 
Fig.|3jc), J will vanish in both small and large U limit. 
The asymptotic behavior goes like J oc Up for Up -C t 
and J oc t 3/2 {U p)^ 1/2 for Up > t. The maximum of J 
is achieved around Up/t — 3 which is of the order 10 _3 i. 
It is still quite challenging to reach such low temperature 
for current cold atom experiment. 

However, in contrast to many systems, at the lowest 
temperature the system is a conventional Bose conden- 
sate which by itself is not the most interesting phase. 
The most interesting phase in this system exists at the 
temperature regime k^T > J. In this regime, the domain 
wall will proliferate and the system will enter a thermal 
mixed state in which all mean-field configurations satis- 
fying condition (i) and (ii) mentioned above are almost 
equally populated. In this thermal mixed state, for each 
site the condensate phase can take 9i = 0, ±27r/3 with 
equal probabilities. Thus, the single boson correlation 
will become short-ranged, i.e. (b\b j) — p s (e~ l6i e l6j ) —> 0. 
However, we notice that as long as the "3-color arrange- 
ment" is satisfied, at every site e 3l9i = 1 is always hold. 
Thus, the triple bosons operator can still possess long- 
range correlation, i.e. (^ 3 6 3 ) = p 3 s {er 3ie ' e 3l8i ) -» p 3 s . 
Since the "3-color arrangement" is enforced by the ki- 
netic energy, thus, it will not be destroyed by thermal 
fluctuations as long as T < t. Therefore we propose that 
within the temperature range J < T < t, there will be a 
"trion superfluid" phase in which the bosons condense in 
triples. This trion superfluid supports many interesting 
properties, such as 1/3 fractionalized vortices, which will 
be studied in the future. 

When temperature is increased to the same order of 
t, the "3-color arrangement" will be destroyed, and the 
long wave-length fluctuation of the U(l) phase will lead 
to a Kosterlitz-Thouless phase transition from the "trion 
superfluid" to the normal state. This transition tempera- 
ture can be estimated from the superfluid stiffness, which 
can be calculated from the free energy response to the 
phase twist. Let b{ = Wie l9i e lq ri , where e lq ri is the ap- 
plied phase twist. In the temperature regime T> J, we 
can ignore the zero-point energy, and the energy func- 
tional for a given {9i} configuration reads E q [9i,Wi] = 

tE {ij) (WiWje- m - ei) e- iq<ri - ri) + h.c.) + f/EiKI 4 , 
and because all the configurations must be taken into ac- 
count under thermal average, the averaged energy func- 
tional reads 

£q[Wi] = ^ E q[6ii W i\ 

= -\ J2(™>je- iq<ri - ri) + h.c) + UJ2 M 4 (4) 

(ij) 1 

It is interesting to note that Eq. Q is the same of a 
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FIG. 4: (Color on line.) (a) Phase diagram. It contains 
three phases: the conventional superfluid phase with boson 
condensation at K point (KSF), three-boson paired (trion) 
superfluid (Trion SF) phase, and the thermal gas phase, (b) 
The long-range correlation for (b\bj) and {SJ 3 S^) as a function 
of temperature ksT/t. (c) Time-of-fiight intensity at T and 
K point as a function of temperature k^T/t. (d-f): Time-of- 
fiight image of momentum distribution for three typical points 
in the phase diagram (a) . Brighter color indicates larger boson 
density. The thin line of honeycomb marks out the Brillouin 
zone. For (b-f), we have set Up/t = 1. 



mean- field energy of bosons in the Kagome lattice without 
frustration, because the sign of hopping is now inverted 
back. Mathematically, it is because taking thermal av- 
erage leads to (e*( ei_e j)) = —1/2. Physically, it can be 
elaborated more clearly by a fractionalization representa- 
tion (35] . Since the kinetic energy frustration is released, 
trion superfluid will have a finite stiffness. It can be ob- 
tained by minimizing with respect to Wi and ex- 
pand the optimal energy in terms of small q, which leads 
to mm{ w .} £ q [wi] ~ — 2t + ^ 4- \tq 2 . Therefore, the 
stiffness is t, which gives an estimate of the Kosterlitz- 
Thouless transition temperature as Tkt — t/ir. 

With the analysis above, we reach the phase diagram 
as shown in Fig.[4|a). In Fig.[4|b), we show the one- 
boson and three-boson correlation function as a function 
of temperature for a fixed Up/t. The one-particle corre- 
lation (b\b 3 ) is calculated as p s £ m e -m-0i) e -m]/T 
Because we expect this correlation vanishes at very low 
temperature of the order J, we can restrict the summa- 
tion to the configurations that satisfy the "3-color ar- 
rangement" since other configuration cost energy of the 
order t and their contribution is negligible in low temper- 
ature. A[9i] is the zero point energy of each condensate. 



While on the other hand, when we calculate the three- 
boson correlation (St 3 S|) = p 3 s e -^-e 3 ) e -E[e t ]/T ^ 

the summation goes over all configurations {9i}, and the 
energy is given by E[6i] = 2t J2uj) cos (^» — fy); where we 
have ignored the zero-point energy since A[9i] is negligi- 
ble compared to E[6j\. Fig.Qb) shows that there indeed 
exists a large temperature window there the one-boson 
correlation function vanishes while the thrcc-boson cor- 
relation function remains finite. 

Detection. Finally we show that the difference between 
these three phases can be detected by a straightforward 
measurement of momentum distribution via TOF image. 
When bosons condense into K points of the Brillouin 
zone, TOF image will display sharp Bragg peaks at K 
point (or its equivalent K' point) and its reciprocal lattice 
vector, as shown in Fig.[4^d) . Moreover, due to the inter- 
ference effect from the phase structure satisfying the "3- 
color-arrangement" , it is easy to show that the strongest 
peaks do not appear in the first Brillouin zone, but in- 
stead at its reciprocal lattice vector points in the second 
Brillouin zone. As temperature increases to the trion 
superfluid phase, all the 3-color arrangements are ther- 
mally mixed, and the Bragg peak at K point disappears 
as the single-boson superfluid order vanishes. However, 
in contrast to normal state, the TOF image is not feature- 
less. In fact, a large honeycomb structure appears and 
the hexagon is 4 times of the area of the first Brillouin 
zone, as shown in Fig.Qe). The intensity at V point is al- 
ways zero, because for zero-momentum component of the 
Fourier transformation is exactly cancelled out due to the 
"3-color arrangement" . At the highest temperature nor- 
mal state, the TOF image becomes a featureless Gaus- 
sian and the intensity at T point becomes the maximum, 
as shown in Fig.[4^f ) . Thus, we predict that the inten- 
sity at K point rapidly decreases at the transition from 
conventional iC-point condensate to trion superfluid, and 
the intensity of T point rapidly increases at the transi- 
tion from trion superfluid to normal state, as shown in 
Fig.^c). Besides, the trion superfluid order can be also 
probed by the 3-point correlation function of the TOF 
image. 
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